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Abstract 

Our aim in this work is to answer, using simulated narrow-band photometry data, the 
following general question: What can we learn about galaxies from these new generation 
cosmological surveys? For instance, can we estimate stellar age and metallicity distribu- 
tions? Can we separate star-forming galaxies from AGN? Can we measure emission lines, 
nebular abundances and extinction? With what precision? 

To accomplish this, we selected a sample of about 300k galaxies with good S/N from the 
SDSS and divided them in two groups: 200k objects and a template library of 100k. We 
corrected the spectra to z — and converted them to filter fluxes. Using a statistical 
approach, we calculated a Probability Distribution Function (PDF) for each property of each 
object and the library. Since we have the properties of all the data from the STARLIGHT- 
SDSS database, we could compare them with the results obtained from summaries of the 
PDF (mean, median, etc). 

Our results shows that we retrieve the weighted average of the log of the galaxy age with 
a good error margin (a w 0.1 — 0.2 dex), and similarly for the physical properties such as 
mass-to-light ratio, mean stellar metallicity, etc. Furthermore, our main result is that we 
can derive emission line intensities and ratios with similar precision. This makes this 
method unique in comparison to the other methods on the market to analyze photometry 
data and shows that, from the point of view of galaxy studies, future photometric surveys 
will be much more useful than anticipated. 

1 Introduction 

This paper is essentially motivated due to the J- PAS project http:/ /j-pas.org, which in a near 
future will produce 8000 sq degrees of imaging in 56 filterbands of 135A up to magnitude 
%ab ~ 24. The survey, besides producing a very high-quality photometric redshift catalog, 
will provide highly informative data important to other astronomical communities as galaxy 
evolution. 
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Figure 1: JPAS 56 filtersystem sensitivity curves. The 51 solid lines are the filters used on 
this work assuming mirror and CCD characteristics and an airmass of 1.3. We also plotted 
in solid fill the SDSS u, g, r, i and z filters in the same conditions. 



Here we describe a bayesian method to "boost the spectral resolution" of these kind of 
photometric data. The idea, which will be described in more detail on section [3j is basically 
that if a high-resolution degrated spectra (or a set of them) is similar to an observed j-spectr4II 
then its physical properties (and even emission lines) would be similar. 

To test the efficiency of this method, we used the STARLIGHT-SDSS [2 J http://starlight.ufsc.br/ 
database as a sandbox. We downloaded the 926246 galaxies where there are measured phys- 
ical properties and emission lines. 



2 From SDSS to JPAS 

Given a SDSS observed spectra, one can convert its energy distribution 0\ to an arbitrary 
observed photometric filter I defined by its transmission curve 2} \ using the simple con- 
version J i = ^ r^r~Hr~ and the error on the filter I, a 2 (Ji), with the relation cr 2 (Ji) = 
J Ti \d\ 

Ai(o"(nA) 2 if^)AA where A; = N\jA\ is the effective filter size and n\ is the spectral error in 
each point. 

The filtersystem curves to JPAS considering an airmass of 1.3 and the expected CCD 
plus Telescope efficiencies are shown on figure [TJ There will be 56 filters, but due to the 
spectral coverage of SDSS we removed the first one in the blue and the four last filters giving 
us an filtersystem of 51 filters plotted in solid lines. For comparsion, we plotted right below 
our filtersystem the SDSS filterset u, g, r, i and z sensitivities through the same air mass. 

The main idea of this work is to obtain high precision properties from low resolution 
data by taking a shortcut from what we measure from high resolution spectra. We consider 



We define a j-spectrum by the set of photometric measurements on J-PAS filter system. 
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Figure 2: Two examples of j-spectrum and its 5 best matches based on y 2 - At the top, in 
green, the two observed spectra examples: an early-type on left panel and a star-forming 
galaxy on the the right. In blue, from top to bottom, the first five best matches. Their 
j-spectrum is represented by the connected red dots. Yellow points are flagged as bad pixels 
and were not considered. Stars represent the points were j-spectrum is not available. 

that if an pair object-template j-spectra are similar, then their measured properties will be 
similar as well. On fig [2] we show two examples of matching two objects with their five best 
matches in respect to our template library librarjU. In green we plotted the SDSS object and 
in blue its four best matches accordingly to their x 2 calculated over their j-spectrum (plotted 
as red dots). On the right panel we have an example of a star- forming galaxy and on the left 
an early-type. 

This idea is similar to the adopted by [5] and [6] and others, but with an important 
difference: Here we do not compare an observed spectrum with a set of models but we 
compare it to a set of another observed spectra. This not only permits us to measure the 
standard physical properties such as mean age and stellar masses that can be measured by 
ohter methodos, but it also allows us to measure indirectly emission lines on data which we 
evidently do not have enough resolution. This is the most relevant advantage of this method. 

3 Method 

To derive the properties p listed down on table [TJ we calculated for 100k objects and two sam- 
ples of comparsion their likelihood functions £{ = e~* L 2 Xi >i . Where ■ = V\ — (On — 

2 In this paper all the galaxies which belongs to the comparsion sample we call library template galaxies. 
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Figure 3: Examples of normalized PDF distributions to age, extinction, Ha and [N n]/H<x 
The left and right boxes corresponds to the left and right observed galaxaies shown on fig. 
[2j In blue dashed line we show the distribution of our base of templates, in magenta the 
likelihood distribution (or the posterior) and the horizontal lines represents in solid black the 
16th and 84th percentiles, in red dashed the average and in solid green the value measured 
by starlight directly on the spectrum. 

aj jBi i) 2 wf 4 a with a scaling factor a-i = l ' ] \M which is determined by interacting 

t—il 1,3 1,1-, 3 

the calculation of a and \ 2 until a convergence critery of Aajj < 10 -5 is accomplished. In 
our simulations, this takes no more than four interactions. The fi term adjusts the width of 
the PDF. It can be adjusted to minimize errors, but we will not treat it here. 

The weighting used to the spectra was defined by wf i j = 

This was selected to have an unbiased weight in the amplitude of the spectra. In other words, 
we would like to give the same importance to the parts with high fluxes (e.g. emission lines) 
than that to that regions that have lower ones (e.g. continuum). 

With the likelihood for each object and base, we calculated as the output of our method 
a PDF estimator. In our case, to estimate the output, we used the likelihood-weighted average 

Z^j Pj£-i,3 

Pi — v 7" — Pout.i 

l^j >~i,3 

4 Sample selection 

Our main sample was retrieved from the STARLIGHT-SDSS database based on very wide 
criteria, trying to get all kinds of objects. Firstly, we separated from the database the 
galaxies which were in the SDSS main galaxy sample [7j, then we selected those who do not 
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have any bad pixel in intervals of 31 A centered on the emission lines Ha, H/3, [N n]A6584, 
[O ll]A3727 and [O iii]A5007 to assure that when we do not measure an emission line it is 
because it is too weak to measure and not because of any observational error. A last selection 
in redshift (0.01 < z < 0.11) was applied in order to have spectral coverage on the wavelength 
interval where JPAS will observe. Those criteria reduced the total number of galaxies on our 
sample from 926246 to 299253 galaxies. All these galaxies were corrected to the rest-frame. 
As a quick-look test, we will not treat the redshift as a variable here and all results will be 
shown to z = 0. 

We then converted all the observed spectra to J-spectra. In case of problems on the 
SDSS spectra (e.g. bad pixels), we changed the observed flux to the best fit from the 
starlight in the cases where the filters have less bad pixels than 50% of the filter width. 
Otherwise, the filter is flagged to be neglected. 

Then, we divided the sample in two sub-samples: one with the galaxies with S/N > 20 
(113821 galaxies) which we call mother library and other one with galaxies with S/N < 20 
(185432 galaxies) which we call object sample. All comparsions made in this paper will be 
in respect of a set of objects and their PDFs calculated over a given library which is a set 
of galaxies from the mother library. From the mother library, we selected two samples of 
galaxies based on two independent diagrams. 

The first one, which we call CMD, is based on the physical analog to the color-magnitude 
diagram to galaxies, the log M* - (log t) diagram. We chopped this diagram in boxes of 0.1 
dex, and on each box we got 10% of the objects distributed along the Ay axis. This gives us 
11952 galaxies on this library. 

The second library, which we call WHAN, is based on the WHAN diagnostic diagram 
introduced by [3] and divide the galaxies basically between star-forming, active nuclei and 
passive or retired galaxies. We did the same cut on this diagram as we did on CMD library, 
but here we changed the Ay to the emission line ratio log[Nn]/Halpha. This gives us 27537 
galaxies on this library. 

5 Results 

To estimate our method's precision, we compared the properties derived by our method 
(output) on the low-resolution j-spectra with the values derived using starlight on the 
high-resolution SDSS spectra (input). So, to do this comparsion, we evaluate, for each 
property p and object i, the Api = ^output -Jfyinput- 

This experiment shows the potential of the proposed method and accomplish with our 
initial objective which is to test the precision of galaxy properties with it. From table [H 
we can resume the precision of our method: It can measure physical properties like age and 
extinction with typical precision of 0.2 dex (or mag, in the case of Ay) and emission lines 
with 0.3 dex and our library selection does not affected the final result, probably because we 
have a number of templates which are in both libraries and/or beacause we have oversampled 
libraries with iVg alaxies of about 10k galaxies. 

Once more, our main result is that we can measure emission lines without having 
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sufficient spectral resolution to do it directly on our data. 



Table 1: Results 



P T i~l in V t" T T 


An 
LAP 


n(ArA (CMT)) 


An 


n(An) fWHAlNH 

U 1 Z_l JJ I I V V M J 


A 

Ay 


0.023 


U.lUb 


0.023 


U.1U1 


(log 


-0.018 


0.199 


-0.013 


n i no 

0.192 


\iog Z,*)l 


-U.U/l 


U. 144 


-u.uzu 


n i a i 

U. 141 


log M/L r 


-0.045 


0.114 


-0.040 


0.110 


log W[ 0jj] 


0.051 


0.223 


0.040 


0.218 


log W H/3 


0.024 


0.145 


0.030 


0.143 


logW [0lIl] 


0.046 


0.245 


0.029 


0.232 


log % Q 


0.010 


0.160 


0.022 


0.157 


log W m] 


-0.028 


0.159 


-0.010 


0.156 


log[N n ]/H a 


-0.045 


0.141 


-0.044 


0.146 


\og[O in ]/Hp 


0.026 


0.250 


-0.000 


0.238 


logHjHf, 


-0.011 


0.107 


-0.010 


0.114 


log S n /H a 


-0.006 


0.172 


-0.016 


0.174 


]og[O n ]/Hp 


0.036 


0.202 


0.016 


0.198 


\og[O in ]/[N n ] 


0.075 


0.265 


0.044 


0.252 



Acknowledgments 

We thank financial support from CNPq, FAPESP, Instituto dc Astronomi'a de Andalucia and the 
CNPq's Instituto Nacional de Ciencia e Tecnologia - Astrofisica. NVA has been supported by CAPES 
(proc. no. 6382-10-0). WS has been supported by Spanish Ministerio de Economi'a y Competitivi- 
dad (grant AYA2010-22111-C03-01) and Brazilian INCT-A. The SEAGal Team wishes to thank all 
researchers involved in the Sloan Digital Sky Survey for their dedication to a project which has made 
the present work possible. 

References 

[1] Bcmtcz, N., 2000, ApJ, 536, 571 

[2] Cid Fernandes R. et al, 2005, MNRAS 358, 363 

[3] Cid Fernandes R. ct al, 2010, MNRAS, 403, 1036 

[4] Driver, S. P, 2011, arXiv:1112.6244 

[5] Gallazzi, A. et al, 2005, MNRAS 362, 41 

[6] Kauffmann, G. et al, D., 2003, MNRAS, 341, 33 

[7] York, D. G et al, 2000, AJ, 120, 1579 



